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PRACTICAL  METHODS  FOR  THE  COMPENSATION  AND 


CONTROL  OF  MULTIVARIABLE  SYSTEMS 


Research  Objectives: 

The  primary  purpose  of  this  research  effort  involved  the  development  of  new  and 
innovative  techniques  for  the  compensation  and  control  of  complex  multivariable  sys¬ 
tems,  such  as  high  performance  aircraft,  ballistic  missiles,  autonomous  vehicles,  and 
industrial  robots.  This  effort  did  not  focus  on  any  one  particular  approach  due  to  the 
diversity  of  assumptions  made  regarding  the  mathematical  and/or  physical  description 
of  the  class  of  systems  considered  as  well  as  the  type  of  performance  required  within 
each  class. 

Fundamental  to  all  of  these  investigations  was  the  practicality  of  the  design  from 
the  point  of  view  of  physical  implementation.  Indeed,  the  revolution  in  computer 
hardware  has  created  an  exceptional  opportunity  for  the  implementation  of  complex 
designs,  especially  through  the  use  of  microprocessors.  Throughout  this  research  effort, 
the  interplay  between  the  theoretical  and  practical  levels  was  of  considerable  impor¬ 
tance. 

To  be  more  specific,  the  simplest  mathematical  model  employed  was  a  linear,  time 
invariant  one  defined  in  either  state  space,  transfer  matrix,  or  polynomial  matrix  terms. 
Greater  emphasis  was  placed  on  polynomial  matrix  models  as  several  recent  investiga¬ 
tions  have  shown  that  such  models  encompass  some  of  the  best  features  of  the  better 

known  transfer  matrix  and  state  apace  m^fc/d^^^^rifjy^o^utational 
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algorithms,  such  as  adaptive  and  block  processing  pole  placement  schemes,  have  been 
developed  for  systems  whose  dynamical  performance  is  described  via  polynomial 
matrices.  Research  efforts  were  directed  at  expanding  theoretical  results  in  the  areas  of 
block  processing  and  adaptive  control  as  well  as  implementing  the  algorithms  on 
devoted,  microprocessor  based  machines. 

The  most  recent  efforts  were  focused  on  the  development  of  multipurpose  controll¬ 
ers  for  a  rather  broad  class  of  systems  whose  desired  performance  can  be  described  in 
“real  world”  Cartesian  space,  but  whose  control  is  achieved  in  an  alternative,  physically 
based  space.  Specific  investigations  in  this  area  were  directed  at  the  control  of  indus¬ 
trial  type  robots.  In  particular,  dynamical  system  simulations  of  various  anthropo¬ 
morphic  arms  were  used  to  develop  a  "dual  drive”  controller  which  enables  a  robot 
equipped  with  force  sensors  to  move  in  any  direction  along  a  smooth  surface  with  a 
specified  force  and  velocity,  irrespective  of  changes  encountered  in  the  surface  contour. 
Here,  the  task  is  defined  in  a  Cartesian  space  while  control  is  implemented  via  joint  axis 
control;  i.e.  in  “joint  space”.  Initial  investigations  in  this  area  have  led  to  the  develop¬ 
ment  of  a  new  computational  solution  to  the  inverse  kinematic  problem  which  could 
have  significant  impact  on  the  practical  implementation  of  advanced  control  designs, 
such  as  the  dual  drive  controller. 

Specific  Research  Accomplishments: 

A.  Multivariable  Adaptive  Control 

The  polynomial  matrix  approach  was  shown  to  represent  a  practical  method  of 
characterizing  the  dynamical  behavior  of  deterministic  systems  from  the  point  of  view  of 
parameter  estimation  for  adaptive  control.  In  particular,  [1]  contains  a  unified 
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could  be  significant  from  the  point  of  view  of  microprocessor  implementation. 

C.  Robotic  Control 

Generally  speaking,  most  industrial  robots  today  possess  a  kinematic  structure 
characterized  by  one  or  more  revolute  joints.  However,  many  of  the  tasks  they  are 
required  to  perform  involve  positioning  in  a  “rectangular”,  Cartesian  space,  such  as 
straight  line  trajectory  motion  between  two  Cartesian  points  or  end  effector  alignment 
with  respect  to  a  defined  point  on  some  planar  surface.  In  such  cases,  control  is  actually 
implemented  at  the  joint  level;  i.e.  in  “joint  space”,  which  will  be  denoted  by  the  n 
(>2)  -dimensional  joint  space  configuration  vector,  §_  ,  while  motion  is  specified  in 
“Cartesian  space”,  which  will  be  denoted  by  the  n-dimensional  Cartesian  space 
configuration  vector  In  view  of  these  observations,  it  is  not  difficult  to  see  that  the 
transformation  problem  from  one  space  to  the  other  represents  perhaps  the  most  funda¬ 
mental  one  in  robotic  manipulation  and  control.  Furthermore,  it  is  of  interest  to  note 
that  this  transformation  question  does  not,  as  yet,  have  a  complete  and  straightforward 
solution  in  the  general  case. 

To  be  more  specific,  the  so-called  forward  kinematic  problem  has  a  known,  general 
solution;  i.e.  it  is  always  possible  to  express  *  in  terms  of  £  in  a  relatively  straightfor¬ 
ward  manner  for  virtually  any  known  kinematic  configuration  through  the  use  of  homo¬ 
geneous  matrix  transformations;  i.e.  a  nonlinear,  n-dimensional  vector  valued  function, 
0.(0)  can  be  obtained  with  relative  ease  such  that 

However,  the  converse  or  inverse  kinematic  problem,  namely  solving  for  a  particular 


discussion  of  direct  and  indirect  strategies  for  parameterizing  multivariable  adaptive 
controllers.  By  considering  unknown,  but  linear  and  time-invariant  systems  in  a  deter¬ 
ministic  setting,  virtually  all  commonly  employed  adaptive  control  strategies  were 
derived  using  pole  placement  notions.  The  particular  problem  which  was  considered  is 
perhaps  the  most  important  parameterization  issue  limiting  the  practicality  of  multivari¬ 
able  adaptive  control,  namely  the  excessive  size  of  the  estimation  problem.  The  polyno¬ 
mial  matrix  formulation  was  shown  to  offer  considerable  advantages  over  other  model¬ 
ling  techniques  from  the  point  of  view  of  reducing  the  number  of  parameters  which  need 
to  be  estimated. 

B.  Block  Processing 

The  general  problem  of  developing  simple  control  algorithms  for  microprocessor 
implementation  motivated  a  closer  look  at  the  effect  which  discrete  periodic  output  con¬ 
trol  laws,  such  as 

ui  =  .  »>1, 

have  on  the  closed  loop  response  of  scalar,  n-dimensional  systems  which  can  be  modelled 
via  the  discrete  equations: 

x1+,  =  Ax-,  -I-  buj ,  y-,  =  ex,. 

It  should  be  noted  that  the  relative  ease  of  implementation  makes  such  control  laws  an 
attractive  alternative  to  conventional  dynamic  compensation  for  process  regulation.  In 
[2],  new  conditions  were  presented  for  the  determination  of  a  sequence  f;  which  produces 
a  dead-beat  response;  i.e.  a  system  output  which  exactly  matches  the  system  input  at 
the  defined  discrete  intervals  after  some  finite  number  of  discrete  steps.  These  results 
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joint  configuration  9  in  terms  of  a  given  Cartesian  configuration  &  is  not  nearly  as 
straightforward. 

In  view  of  the  preceding,  the  main  purpose  of  [3]  was  to  present  a  new,  computa¬ 
tional  (numerical)  solution  to  the  “general  version”  of  the  inverse  kinematic  problem.  In 
particular,  it  was  shown  in  (3)  that  a  relatively  simple  dynamical  system  can  be  con¬ 
structed  which,  when  “driven”  by  a  desired,  time-varying,  Cartesian  configuration  vec¬ 
tor,  2Ld(t)  ,  produces  not  only  the  corresponding  fl^t)  ,  such  that  £L(i?d(t))  =  ^(t)  ,  but 
also  ^j(t)  and,  if  required,  i^t)  ,  as  well.  It  should  be  noted  that  the  ability  to  solve  the 
general  inverse  kinematic  problem  in  real  time  is  fundamental  to  the  implementation  of 
certain  more  advanced  control  algorithms,  such  as  the  dual  drive  controller  noted  ear¬ 
lier. 

D.  Microprocessor  Implementation 

As  stated  at  the  outset,  one  of  the  primary  objectives  of  our  overall  research  effort 
was  to  develop  control  algorithms  which  can  be  physically  implemented  in  a  rather 
direct  manner  via  microprocessors.  In  light  of  this  objective,  research  within  the 
Laboratory  for  Engineering  Man/Machine  Systems  (LEMS)  at  Brown  University  focused 
on  the  development  of  a  general  purpose,  68000  based  microprocessor  controller  which 
could  be  employed  in  a  variety  of  control  environments.  At  present,  a  preliminary  ver¬ 
sion  of  such  a  unit  has  been  constructed  and  is  being  used  to  control  a  single  axis  of  an 
IBM  RS/1  Cartesian  robot.  This  research  is  directed  at  not  only  testing  the  micropro¬ 
cessor  controller,  but  also  at  developing  and  testing  certain  of  the  control  algorithms 
which  have  evolved  under  AFOSR  sponsorship.  Reference  [4]  will  contain  a  detailed 


could  be  significant  from  the  point  of  view  of  microprocessor  implementation. 

C.  Robotic  Control 

Generally  speaking,  most  industrial  robots  today  possess  a  kinematic  structure 
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However,  the  converse  or  inverse  kinematic  problem,  namely  solving  for  a  particular 


discussion  of  this  particular  part  of  our  general  research  effort. 
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joint  configuration  9  in  terms  of  a  given  Cartesian  configuration  &  is  not  nearly  as 
straightforward. 

In  view  of  the  preceding,  the  main  purpose  of  [3]  was  to  present  a  new,  computa¬ 
tional  (numerical)  solution  to  the  “general  version”  of  the  inverse  kinematic  problem.  In 
particular,  it  was  shown  in  (3)  that  a  relatively  simple  dynamical  system  can  be  con¬ 
structed  which,  when  “driven”  by  a  desired,  time-varying,  Cartesian  configuration  vec¬ 
tor,  2Ld(t)  ,  produces  not  only  the  corresponding  fl^t)  ,  such  that  £L(i?d(t))  =  ^(t)  ,  but 
also  ^j(t)  and,  if  required,  i^t)  ,  as  well.  It  should  be  noted  that  the  ability  to  solve  the 
general  inverse  kinematic  problem  in  real  time  is  fundamental  to  the  implementation  of 
certain  more  advanced  control  algorithms,  such  as  the  dual  drive  controller  noted  ear¬ 
lier. 

D.  Microprocessor  Implementation 

As  stated  at  the  outset,  one  of  the  primary  objectives  of  our  overall  research  effort 
was  to  develop  control  algorithms  which  can  be  physically  implemented  in  a  rather 
direct  manner  via  microprocessors.  In  light  of  this  objective,  research  within  the 
Laboratory  for  Engineering  Man/Machine  Systems  (LEMS)  at  Brown  University  focused 
on  the  development  of  a  general  purpose,  68000  based  microprocessor  controller  which 
could  be  employed  in  a  variety  of  control  environments.  At  present,  a  preliminary  ver¬ 
sion  of  such  a  unit  has  been  constructed  and  is  being  used  to  control  a  single  axis  of  an 
IBM  RS/1  Cartesian  robot.  This  research  is  directed  at  not  only  testing  the  micropro¬ 
cessor  controller,  but  also  at  developing  and  testing  certain  of  the  control  algorithms 
which  have  evolved  under  AFOSR  sponsorship.  Reference  [4]  will  contain  a  detailed 
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